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ABSTRACT 

The experimental emission probabilities of complex fragments by low energy com- 
pound nuclei and their dependence upon energy and Z value are compared to the 
transition state rates. Intermediate-mass- fragment multiplicity distributions for a va- 
riety of reactions at intermediate energies are shown to be binomial and thus reducible 
at all measured transverse energies. From these distributions a single binary event 
probability can be extracted which has a thermal dependence. A strong thermal sig- 
nature is also found in the charge distributions. The n-fold charge distributions are 
reducible to the 1-fold charge distributions through a simple scaling dictated by fold 
number and charge conservation. 



1. Transition State Rates and Complex Fragment Decay Widths 

The rates for fission decay, as well as for chemical reactions, are calculated most 
often by means of the transition state method . In this approach, the reaction rate 
is equated to the flux of phase space density across a "suitably" located hyperplane 
normal to the "reaction coordinate". The "suitable" location is typically chosen at 
a saddle point in collective coordinate space, which corresponds to a bottleneck in 
phase space. A smart choice of the transition state location should minimize the 
number of phase space trajectories doubling back across the hyperplane. 

The surprising success of the transition state method has prompted attempts 
to justify its validity in a more fundamental way, and to identify regimes in which 
deviations might be expected |]^, ^ Q . In what follows we shall compare experimental 
decay rates for complex fragment emission with transition state predictions, and 
search for energy E and atomic number Z dependent deviations that can be expected 
to exist. 

The transition state expression for the flssion decay width is: 



(1) 



where p{E) is the level density of the compound nucleus, p*{E — Bf — e) is the level 
density at the saddle point, Bf is the fission barrier, e is the kinetic energy over the 
saddle along the fission coordinate and 



1 _ d[\np*{x)] 
Yf ~ dx 



(2) 



For the one dimensional case in which the only degree of freedom treated explicitly 
is the reaction coordinate, the decay width takes the form: 



where T is the temperature of the transition state. Now both level densities corre- 
spond to the same number of degrees of freedom. The quantity huj is the oscillator 
phonon associated with the ground state minimum. 

The emission of complex fragments can be treated in an analogous fashion by 
introducing the ridge line of conditional saddle points [^. Each mass or charge 
emission can be associated with a conditional barrier. These barriers can be measured 
with techniques similar to those used to determine fission barriers . Recently, nearly 
complete ridge lines have been determined for several nuclei: ''^Br [0, ^°'^^Mo and 

The emission rate of a fragment of a given mass or charge can still be described 
by an expression similar to that of Eq. @. The quantity Bf becomes the conditional 
barrier Bz] but what is now the meaning of hujl Is there a single value of hu for all 
the channels or has each channel its own characteristic frequency? We shall endeavor 
to answer this question experimentally. 

An additional aspect of the problem has been studied by Kramers in his seminal 
work 0. Kramers considered the diffusion of the system from the reactants' region to 
the products' region from the point of view of the Fokker-Planck equation. The new 
parameter entering the problem is the viscosity coefficient, which couples the reaction 
coordinate to the heat bath. The stationary current solution found by Kramers leads 
to expressions for the reaction rates similar to that of the transition state theory, 
differing only in the pre-exponential factor, which now includes the viscosity. More 
recent work on the same equation has shown that if the system is forced to start at 
time t=0 at the ground state minimum, a transient time Tf exists during which the 
reaction rate goes from zero to its stationary value 0. Both effects would decrease 
the overall fission rate compared to the transition state prediction. 

These effects have been advocated as an explanation for the large number of 
prescission neutrons observed in the fission of many systems p|, |10|, |ll], |12|, 0, in 
apparent contradiction with the predictions of the transition state method |TT|, 112 . 
The prescission neutrons can be emitted either before the system reaches the saddle 
point, or during the descent from saddle to scission. Only the former component. 
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however, has any bearing on possible deviations of the fission rate from its transition 
state value, and the separation of the two components is very difficult indeed. 

Furthermore, it has been suggested that the viscosity and the transient time may 
depend on the collectivity of the reaction coordinate More specifically, the 



reaction coordinate for a very asymmetric decay should have little collectivity, while 
that for a symmetric decay should be very collective. A study of prescission particles 
as a function of the size of the emitted fragment claims to have observed such and 



effect 11, 14, 15 



We are going to show that the presence or absence of the effects discussed above 
can be directly observable in the excitation functions for the emission of fragments 
with different Z values. Our procedure uses the transition state prediction as a null 
hypotheses, and involves only the replotting of experimental data without using any 
specific model. The cross section for the emission of a fragment of a given Z value 
can be written as: 

Oz = fTo— = (To , (4) 

J-T in + J-pH---- 

where ctq is the compound nucleus formation cross section and Tt, r„, Tp, V z are 
the total-, neutron-, proton-, and Z-decay widths, respectively. Notice that Vt is 
essentially independent of Z if we confine our observations to the excitation energy 
region where the complex fragment emission probability is small. 
We now rewrite Eq. as follows: 

a^^^2ME_m . _ _ E-.), (5) 

where Tz is the temperature at the conditional saddle point and E-l.^, ^ the 
ground state and saddle point rotational energies. In this way, the left hand side of 
the equation contains the complex fragment cross section which can be measured, 
and other calculable quantities that do not depend on Z, except Tz which is only 
weakly dependent on Z. The right hand side contains only the level density at the 
conditional saddle calculated at the intrinsic excitation energy over the conditional 
saddle, which is calculable if the barrier height is known. 

By using the standard Fermi gas level density expression, one can rewrite Eq. (§) 
in the following way which takes out the A-dependence of the level density: 



In 



■f 




-Bz-E^), (6) 

where az, ctn are the saddle and ground state level density parameters. A plot of the 
left hand side of this equation versus the square root of the intrinsic excitation energy 
over the saddle should give a straight line, and the slope of this straight line should 
give the square root of a^/a^. 

Recently, the excitation functions for a large number of fragment Z values have 
been measured for the following systems: ^^Br 0, ^°'^^Mo 0, and ^^^^^^"^In 0. The 
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Fig. 1: a) The excitation functions (cross sections vs excitation energy) for complex fragments of 
some typical Z values emitted from the compound nucleus ^^Mo produced in the reaction ^^Kr+^^C 
at beam energies ranging from 6.2 to 12.2 MeV/u. b) The az/a,n values and c) conditional barriers 
Bz, both extracted by fitting the excitation functions with a transition state formalism. The solid 
lines in a) correspond to the fit using an energy level parameter a„ = A/8. 
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Fig. 2: \n Rf /2y/a^ (see Eq. and text) vs the square root of the intrinsic excitation energy for 
four compound nuclei: ^^Br a), ^"Mo b), ^"^Mo c), and ^^'^'^^'^In d). All the excitation functions for 
the indicated Z range are included for each compound nucleus. The solid lines are the linear fits to 
the data. The error bars are smaller than the size of symbols. 

corresponding conditional barriers have been extracted by fitting the excitation func- 
tions with the transition state formahsm. A level density parameter = A/ 8 was 
assumed in the fitting. As an example, Fig. |l]a shows some of the excitation functions 
for fragments with Z-values from 5 to 25 for the compound nucleus ^^Mo. The solid 
lines in Fig. ^ correspond to the best fit to the experimental data. The extracted ra- 
tios az/cLn are close to unity for all Z values (see Fig. Pd). The extracted conditional 
barriers increase from 30-45 MeV as the charge of the emitted fragment increases (see 
Fig. Ic). 

Eq. (H) suggests that it should be possible to reduce all the excitation functions for 
the emission of different complex fragments from a given system to a single straight 
line. In Fig. ^ all the excitation functions associated with each of four compound 
nuclei (''^Br, ^°Mo, ^^Mo, and '^^^'^'^'^In) are plotted according to Eq. (^). There are 
20, 21, 21, and 9 excitation functions for ''^Br, ^°Mo, ^^Mo, and ^^°'^^^In, respectively. 
We see that all the excitation functions for each Z- value fall with remarkable precision 
on the same line that is in fact straight, has a slope near unity and passes closely 
through zero. 

As a final virtuoso touch, we can try to collapse ALL the excitation functions for 
ALL Z values and for all compound nuclei into a single straight line. The resulting 
plot for four different compound nuclei is shown in Fig. 0. It includes a total of 
71 excitation functions spanning a Z range from 3 to 25. The collapse of all the 
experimental excitation functions for all of the different Z-values and all the systems 
onto a single straight line is strong evidence for the validity of the transition state 
formalism and for the absence of Z- and ii^-dependent deviations. In particular, one is 
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led to the following conclusions: a) Once one removes the phase space associated with 
the non reactive degrees of freedom at the conditional saddle point, the reduced rates 
are IDENTICAL for fragments of all Z-values. Within the experimental sensitivity, 
the quantity huj in Eq. @ appears to be Z independent, b) For all fragments, there 
is no deviation from the expected linear dependence over the excitation energy range 
from 50-130 MeV explored. This seems to rule out, for all Z- values, transient time 
effects which should become noticeable with increasing excitation energy, c) The 

slope, which corresponds to the y^a^/a^, is essentially 1 for all Z values of all systems 
studied, d) The intercept of the straight line, which is associated with the channel 
frequency tu, is essentially zero and shows no obvious dependence on the fragment 
Z- values (i.e., the collectivity). 

We conclude that in this extended data set there is no evidence for transient effects 
either directly or through their expected dependence upon the mass of the emitted 
fragment. Furthermore it appears that the channel frequency is the same for all the 
different Z decay channels. 



2. Reducibility at intermediate energies 

At low excitation energies, complex fragments are emitted with low probability 
by a compound nucleus mechanism ^ . At increasingly larger energies, the prob- 
ability of complex fragment emission increases dramatically, until several fragments 



are observed within a single event |T8| , |T9|, ^ . The nature of this multifragmentation 
process is at the center of much current attention. For example, the timescale of 
fragment emission and the associated issue of sequentiality versus simultaneity are 
the objects of intense theoretical and experimental study. 
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Recent experimental work ||21|, g^] has shown that the excitation functions for 
the production of two, three, four, etc. fragments give a characteristically linear 
Arrhenius plot, suggesting a statistical energy dependence. 

A fundamental issue, connected in part to those mentioned above, is that of re- 
ducibility: can multifragmentation be reduced to a combination of (nearly) indepen- 
dent emissions of fragments? More to the point, can the probability for the emission 
of n fragments be reduced to the emission probability of just one fragment? 

Recently, it has been experimentally observed in many reactions that for any value 
of the transverse energy Et, the n-fragment emission probability P„ is reducible to 



the one-fragment emission probability p through a binomial distribution ^ 



777 I 

nl[m — n)l 

This empirical evidence indicates that multifragmentation can be thought of as a 
special combination of nearly independent fragment emissions. The binomial combi- 
nation of the elementary probabilities points to a combinatorial structure associated 
with a time-like or space-like one-dimensional sequence. It was also found that the 
log of such one-fragment emission probabilities (logp) plotted vs l/y/E^ (Arrhenius 
plot) gives a remarkably straight line. This linear dependence is strongly suggestive 
of a thermal nature for p, 

p = e-«/^ (8) 

under the assumption that the temperature T oc ^/E* where E* is the excitation 
energy. Examples of the binomial decomposition of the ra-fragment emission proba- 
bilities Pn into a one-fragment emission probability p, and the resulting Arrhenius plot 
for p is given in Fig. ^. The extraordinary quantitative agreement between the calcu- 
lations and the experimental data confirms the binomiality of the multifragmentation 
process. 

The more directly interpretable physical parameter contained in this analysis is 
the binary barrier B (proportional to the slope of the Arrhenius plot in Fig. |^). One 
may wonder why a single binary barrier suffices, since mass asymmetries with many 
different barriers may be present. Let us consider a barrier distribution as a function 
of mass asymmetry x of the form B = Bo + ax^, where Bo is the lowest barrier in the 
range considered. Then, 



nujQ J \a J 



p 



Thus the simple form of Eq. (||) is retained with a small and renormalizable pre- 
exponential modification. 

One possible interpretation of the reducibility discussed above is sequential decay 
with constant probability p. Assuming that the (small) fragments, once produced 



7 



E/A = 50 MeV 




Fig. 4: For the -^^^Xe induced reactions at E/A=50 MeV (left figure) and the '^^Ar+^^''Au reactions 
at E/A=50, 80 and 110 MeV (right figure): (left panel) the reciprocal of the single fragment emission 
probability 1/p as a function of 1/ a/EJ"; and (right panel) the parameter m (number of the throws 
in the binomial distribution) and the probability P(n) of emitting n intermediate mass fragments 
(IMF, 3 < Z < 20) as a function of the transverse energy Et. The solid lines through the excitation 
functions correspond to binomial distributions calculated with the given values of m and p. (See 
Eq. ®) 
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do not generate additional fragments or disappear, the binomial distribution follows 
directly. In this framework, it is possible to translate the probability p into the mean 
time separation between fragments. In other words, we can relate the n-fragment 
emission probabilities to the mean time separations between fragments. The validity 
of this interpretation is testable by experiment. 

We have tried to find alternative explanations to the sequential description for 
the binomial distributions with thermal probabilities. An obvious model is a chain 
of m links with probability p that any of the links is broken. The probability that n 
links are broken is given by Eq. (|^). This result is, of course, strictly dependent on 
the dimensionality of the model, and its relevance to multifragmentation is unclear. 
Nevertheless, it stresses again the fundamental reducibility of the multifragmentation 
probability to a binary breakup probability p. 



3. Charge Distributions 

These aspects of reducibility and thermal scaling in the integrated fragment emis- 
sion probabilities lead naturally to the question: Is the charge distribution itself re- 
ducible and scalable? In particular, what is the charge distribution form that satisfies 
the condition of reducibility and of thermal dependence? 

Let us first consider the aspect of reducibility as it applies to the charge distri- 
butions. In its broadest form, reducibility demands that the probability p{Z), from 
which an event of n fragments is generated by m trials, is the same at every step 
of extraction. The consequence of this extreme reducibility is straightforward: the 
charge distribution for the one-fold events is the same as that for the n-fold events 
and equal to the singles distributions, i.e.: 



P(1)(Z) = P(„)(Z) = Ps^ngles{Z) = p{Z) . 



(10) 



We now consider the consequences of the thermal dependence of p p3[ on the 
charge distributions. If the one-fold = ra-fold = singles distributions is thermal, then 



B(Z) 

P{Z) oc — 



or TlnP(Z) oc —B{Z). This suggests that, under the usual assumption Et oc E* 
23| , the function 

^,\tiP{Z) = D{Z) (12) 

should be independent of Et- 

In the '^^Ar+^^''Au reaction, as in other reactions |2^, |2^, the IMF charge distri- 
butions are empirically found to be nearly exponential functions of Z 



Pn{Z) OC e 



(13) 
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"Ar+^^^Au, E/A=110 MeV 
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Fig. 5: The exponential fit parameter a„ (from fits to the charge distributions, see Eq.(p^) is plotted 
as a function of 1/y/Rt. The solid lines are a fit to the values of a„ using Eq. ([l5|). 



In light of the above considerations, we would expect for q;„ the following simple 
dependence 

1 1 

(14) 



a„ oc — oc 



for all folds n. Thus a plot of a„ vs 1/ \/E't should give nearly straight lines. This is 
shown in Fig. | for ^eAr+^^^Au at E/A=llQ MeV. 

The expectation of thermal scaling appears to be met quite satisfactorily. For each 
value of n the exponent a;„ shows the linear dependence on 1/ ^/El anticipated in Eq. 
(0). On the other hand, the extreme reducibility condition demanded by Eq. ([To|) , 
namely that ai = ^2 = ... = a„ = a, is not rigorously met. Rather than collapsing 
on a single straight line, the values of q;„ for the different fragment multiplicities are 
offset one with respect to another by what appears to be a small constant quantity. 

In fact, one can fit all of the data remarkably well, assuming for a„ the form: 

OLn = —F= + nC (15) 



which implies: 



Or, 



/Ei 



K 

— + nc 



(16) 



or more generally, for the Z distribution: 

Pn{Z) oc e-^-"^^. (17) 

Thus, we expect a more general reducibility expression for the charge distribution of 
any form to be: 

[\YiPn{Z)+ncZ]JEt = F{Z) (18) 
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for all values of n and Ef. This equation indicates that it should be possible to reduce 
the charge distributions associated with any intermediate mass fragment multiplicity 
to the charge distribution of the singles. 

What is the origin of the regular offset that separates the curves in Fig. ^? The 
general form of Eq. (|1^) suggests the presence of an entropy term that does not 
depend explicitly on temperature. The general expression for the free energy (in 
terms of enthalpy H, temperature T and entropy S) 

AG = AH{Z) - TAS{Z) (19) 

leads to the distribution 

P{Z) oc e (20) 

Typically, AS* is of topological or combinatorial origin. For instance, a factor of 
this sort would appear in the isomerization of a molecule involving a change of sym- 
metry. In our specific case AS* may point to an asymptotic combinatorial structure 
of the multifragmentation process in the high temperature limit. As an example, we 
consider the Euler problem of an integer to be written as the sum of smaller integers, 
and calculate the resulting integer distribution. Specifically, let us consider an integer 
Zq to be broken into n pieces. Let nz be the number of pieces of size Z. It can be 
shown that 



nz = —e'^ = cn'e-""^. (21) 



This expression has the correct asymptotic structure for T — >■ oo required by Eq. ([17|) . 
The significance of this form is transparent: First, the overall scale for the fragment 
size is set by the total charge Zq. Second, for a specific multiplicity n, the scale is 
reduced by a factor n to the value Zo/n. 

4. Phase Coexistence 



While Eq. ( pTf ) obviously implies charge conservation, it is not necessary that 
charge conservation be implemented as suggested by it. In fact it is easy to envisage 
a regime where the quantity c should be zero. Sequential thermal emission is a case 
in point. Since any fragment does not know how many other fragments will follow 
its emission, its charge distribution can not reflect the requirement of an unbiased 
partition of the total charge among n fragments. We have in mind a liquid drop 
evaporating fragments of different size and binding energy. "Charge conservation" 
will affect the distribution minimally, unless evaporation consumes the entire system, 
and even then, not in the sense of an unbiased partition. 

On the other hand, in a simultaneous emission controlled by a n-fragment transi- 



tion state fragments would be strongly aware of each other, and would reflect 

such an awareness through the charge distribution. 

The question then arises whether c = or c > 0, or even better, whether one 
can identify a transition from a regime for which c = to a new regime for which 
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'Ar + ""Au, E/A=110 MeV 
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Fig. 6: Top panel: the n gated charge distributions Pn{Z) for the reaction "^^Ar+^^^Au at E/A—\\Q 
MeV. The charge distributions were constructed from events with £'t=650±20 MeV and n—l-b. 
Middle panel: the "reduced" charge distribution ||2^ for the same data using the indicated value of 
c. (The data here are normalized at Z=3). Bottom panel: the log of the ratio of P2{Z) / P^{Z) . The 
slope corresponds to c for n^2 (see Eq.(|2^)). The statistical error bars are shown for errors larger 
than the symbol size. 



c > 0. In order to answer this question, we have studied the charge distributions as a 
function of fragment multiphcity n and transverse energy for a number of systems 
and excitation energies. Specifically, we will present data for the reaction ^^Ar+^^^Au 
at E/A=m and 110 MeV and the reaction ^^^Xe+^^^Au at E/A=hQ and 60 MeV. 

A general approach for measuring c, which does not depend on any specific form 
for the charge distribution, is to construct at each Et the ratio 

Pn{Z) z 



A value of c can be extracted for each n by taking the log of this ratio and finding the 
slope of the resulting graph (see bottom panel of Fig. |^). A weighted average (over all 
IMF multiplicities n) for c can then be constructed at all Et. Alternatively, a can 
be constructed in terms of the differences in FiZ^ (see Eq. ([T8|) ) between any pairs of 
n values and minimized as a function of c. These procedures yield essentially the same 
results. These results are reported in Fig. |^ for the ^^^Xe+^^''Au and ^^Ar+^^'^Au 
reactions. 

It is interesting to notice that for all reactions and bombarding energies the quan- 
tity c starts at or near zero, it increases with increasing Et for small Et values, and 
seems to saturate to a constant value at large Et. 

This behavior can be compared to that of a fluid crossing from the region of 
liquid-vapor coexistence to the region of overheated and unsaturated vapor. In the 
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Fig. 7: Plots of the coefRcient c versus Et for the reactions ^^^Xe+^^^Au at i5/A=50 and 60 MeV 
(top panel) and ^^Ar+^^'^Au at E/A^80 and 110 MeV (bottom panel). The error bars are statistical. 



coexistence region, the properties of the saturated vapor cannot depend on the total 
mass of fluid. The presence of the liquid phase guarantees mass conservation at all 
average densities for any given temperature. A change in mean density (volume) 
merely changes the relative amount of the liquid and vapor, without altering the 
properties of the saturated vapor. Hence the vapor properties, and, in particular, the 
cluster size distributions cannot reflect the total mass or even the mean density of 
the system. In our notation, c = 0. 

On the other hand, in the region of unsaturated vapor, there is no liquid to insure 
mass conservation. Thus the vapor itself must take care of this conservation, at least 
grand canonically. In our notation, c > 0. 

This description should not be taken too literally, for a variety of reasons, one of 
which is the finiteness of the system. The c = regime may signify an evaporative-like 
emission from a source which survives as a charge conserving residue (liquid), while 
the c > regime may signify the complete vaporization of the source. 

In order to test these ideas for a finite system, percolation calculations were 
performed for systems of Zo=97, 160 and 400 as a function of the percentage of bonds 
broken (pf,) in the simulation. Values of c were extracted (using Eq.(^)) as a function 
of Pb- 

The results are shown in Fig. |^. Guided by the insight gleaned from the approxi- 
mate solution to Euler's problem (see Eq. (|2T|) ) we have scaled the extracted values of 



13 



1^ 




0.5 0.6 0.7 0.8 0.9 1.0 



Pb 

Fig. 8: A plot of cZo versus the percentage of broken bonds pb from a percolation calculation 
for three systems Zq~%7 (circles), Zq—\QQ (squares) and Zo=400 (diamonds). The statistical error 
bars are shown for errors larger than the symbol size. 



c by the source size Zq in order to remove this leading dependence and to evidentiate 
the true finite size effects. For values of pb smaller than the critical (percolating) 
value (p^"* ~ 0.75 for an infinite system), we find c = 0. This is the region in which 
a large (percolating) cluster is present. As goes above its critical value, the value 
of c increases, and eventually saturates in a way very similar to that observed exper- 
imentally. Due to the finiteness of the system the transition is smooth rather than 
sharp and can be made sharper by increasing the size of the system. 

Before proceeding, let us remind ourselves that charge conservation is not a finite- 
size effect. For instance, the chemical potential, introduced in statistical mechanics 
to conserve mass, survives the thermodynamical limit and retains its meaning for an 
infinite system, despite the fact that the extensive thermodynamic quantities go to 
infinity. In our case, while it is true that c goes to zero or that 1/c goes to infinity, it 
is also true that the product cZq tends to a finite limit nearly independent of Zq. 

The significance of the actual experimental value of c in the region where it seems 
to saturate is unclear. In Eg. (pTD , c takes a direct meaning for the Euler problem: 
c = 1/Zq. It should be noted that our analysis is not directly comparable to the Euler 
solution (Eq.(|2T|)) since we have restricted ourselves to a limited region (3 < Z < 20) 
of the total charge distribution for our study of how the source is partitioned into 
different IMF multiplicities. It must also be appreciated that Eq. (^) and the as- 
sociated dependence of c upon Zq are characteristic of a one-dimensional percolation 
model. In light of the points mentioned above, it is not unexpected that c appears 
to be proportional, but not equal, to 1/Zq in the three-dimensional percolation cal- 
culation reported in Fig. |. An interpretation of c in terms of the source size may 
be possible when more data and a better understanding of the percolation of finite 
systems are available. 
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